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^ ■ Abstract 



We introduce a new parameterization of four-fermion operator matrix elements 
which does not involve quark masses and thus allows a reduction of systematic 
uncertainties. In order to simplify the matching between lattice and continuum 
renormalization schemes, we express our results in terms of renormalization group 
invariant i?-parameters which are renormalization-scheme and scale independent. 
As an application of our proposal, matrix elements of A/ = 3/2 and SUSY AS" = 2 
operators have been computed. The calculations have been performed using the 
tree-level improved Clover lattice action at two different values of the strong cou- 
pling constant (/? = G/g^ = 6.0 and 6.2), in the quenched approximation. Renor- 
malization constants and mixing coefficients of lattice operators have been ob- 
tained non-perturbatively. Using lowest order xPT, we also obtain {7nT\Oj\K)^^^ = 
(0.11 ± 0.02) GeV^ and {TT7r\0s\K)^P^ = (0.51 ± 0.05) GeV^ at ^ = 2 GeV. 
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1 Introduction 

Since tlie original proposals of using lattice QCD to study hadronic weak de- 
cays [l]-[3], substantial theoretical and numerical progress has been made: the 
main theoretical aspects of the renormalization of composite four-fermion op- 
erators are fully understood [4,5]; the calculation of K^-K^ mixing, expressed 
in terms of the so-called renormalization group invariant 5-parameter Bk-, has 
reached a level of accuracy which is unpaired by any other approach [6]- [9]; 
increasing precision has also been gained in the determination of the electro- 
weak penguin amplitudes necessary to the prediction of the CP-violation pa- 
rameter e'/e [10]-[16]; attempts to compute the matrix element of the QCD 
penguin operator Oq exist [17]. Finally, matrix elements of AS* = 2 operators 
which are relevant to study FCNC effects in SUSY models have been also 
computed [16,18]. 

Following the common lore, matrix elements of weak four-fermion operators 
are given in terms of the so-called 5-parameters which measure the deviation 
of their values from those obtained in the Vacuum Saturation Approximation 
(VSA). A classical example is provided by the matrix element of the AS* = 2 
left-left operator O^^^"^ = S7^(l —75)^ S7^(l —75)^ relevant to the prediction 
of the CP-violation parameter e 

{K'\0^'='\K') = iMlflBK . (1) 

VSA values and 5-parameters are also used for matrix elements of AS* = 1 
operators entering e'/e, in particular Oq = Sal/j.il - l5)d/3T,qQf5liJ.{'^ + lb)<la 
and Og = 3/2sc,7/,(l - 75)ci/3 Eg (^qQulixi^ + IblQa [19,20] 
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Contrary to {K^\0^^^'^\K^) and (7r7r|06|i^)/=o which vanish in the chiral 
limit, the matrix element of the left-right operator (7r7r|08|-ft')/=2 remains fi- 
nite. For this reason, the dependence of this amplitude on quark masses is 
expected to be smooth because it only enters in higher orders of the chiral 
expansion 

/ M^ \^ 
-^^ ^Ci + C2{ms + ma) + 0{ml), (3) 



where Ci,2 are constants, independent of quark masses. 

Quark masses, however, appear exphcitly in (2). Since in the VSA (and in 
the 1/A^ expansion [21]) the expression of the matrix elements is quadratic in 
^™s + "fT^d, predictions for the physical amplitudes are heavily affected by the 
specific value taken for this quantity. Contrary to fx, M^, etc., quark masses 
are not directly measured by experiments and the present accuracy in their 
determination is still rather poor [22]. Therefore, the "conventional" param- 
eterization (2) introduces a large systematic uncertainty in the prediction of 
the physical amplitudes of {Oq)j=o and {Os)i=2 (and of any other left-right 
operator). Moreover, whereas for 0^^^^ we introduce B^ as an alias of the 
matrix element, by using (2) we replace each of the matrix elements with 2 
unknown quantities, i.e. the 5-parameter and nis + rrid. Finally, in many phe- 
nomenological analyses, the values of the 5-parameters of {Oq)i=o and {Os)i=2 
and of the quark masses are taken by independent lattice calculations, thus 
increasing the spread of the theoretical predictions Q. All this can be avoided 
in the lattice approach, where matrix elements can be computed from first 
principles. 

In this paper, we propose a new parameterization of matrix elements in terms 
of well known experimental quantities, without any reference to the VSA and 
therefore to the strange (down) quark mass. This results in a determination of 
physical amplitudes with smaller systematic errors. As an application of our 
proposal we have reanalyzed the lattice correlation functions considered in [16] 
to estimate matrix elements of A/ = 3/2 operators and of the operators of the 
most general AS = 2 Hamiltonian. By comparing the results of the present 
study with those of ref. [16] we show all the advantages of the new parame- 
terization. We give the results for operators renormalized non-perturbatively 
in the RI (MOM) scheme [14,16,23,24]. 

We also introduce a Renormalization Group Invariant (RGI) definition of 
matrix elements (and Wilson coefficients). This definition generalizes to an 
arbitrary basis a concept which has been adopted very successfully for the 
K^-K^ mixing amplitude, which is usually written in terms of the RGI B- 
parameter, Bx- With our definition, Wilson coefficients and operators are 
renormalization-scale and scheme independent. This will hopefully avoid some 
confusion existing in the literature. This confusion is generated by the fact 
that (perturbative) Wilson coefficients and (non-perturbative) matrix ele- 
ments have been computed using different techniques, regularizations, renor- 
malization schemes (different versions of NDR and HV, DRED, RI-MOM with 
different external states) and renormalization scales. In this way, the effective 
Hamiltonian is splitted in terms which are individually scheme and scale in- 



^ We will discuss the correlation between the value of the S-paranieter and the 
quark masses in sec. 4. 



dependent. 

The remainder of the paper is organized as follows: in sec. 2 we define operators 
and matrix elements considered in the present study and introduce the new 
parameterization of the matrix elements; in sec. 3 we give the RGI definition 
for a generic operator basis; in sec. 4 we present our new numerical results and 
compare them with those obtained with the "conventional" definition of the B- 
parameters; in sec. 5 we give our best estimates of the matrix elements and in 
sec. 6 we present our conclusions. All details concerning the non-perturbative 
renormalization of lattice operators and the extraction of matrix elements 
from correlation functions are not discussed here since they can be found in 
refs. [14,16]. In particular, in ref. [16], the same set of lattice data was analyzed. 
We compare the results of the present study to those of this reference. 

2 Matrix elements without quark masses 

In this section we introduce the notation and define operators and matrix el- 
ements used in this paper. The new parameterization, that does not involve 
any reference to quark masses, is defined here. An alternative parameteriza- 
tion which (for reasons explained below) has not been used in our numerical 
analysis, but may be useful in the future, is also considered in this section. 

• A^ = 2 operators: 



The analysis of K^ — K^ mixing with the most general AS* = 2 effective 
Hamiltonian requires the knowledge of the matrix elements {K^\Oi\K^) of 
the following operators [25]-[27] 

Oi = s"7^(l - 75)rf" Am(1 - 75)c^^ 

02 = s"(l-75)rf"s''(l-75)rf^ 

03 = s"(l-75)rf^s'^(l-75)rf", (4) 

04 = s-(l-75)rf"s'^(l+75)rf^ 

o^ = r{i-^,)d^-s^{i + ^,)r, 

where a and (3 are color indices and O^^^"^ = Oi. For {K'^\Oi\K'^) , only the 
parity-even parts of the operators of eqs. (4) contribute. {K^\Oi\K'^) is ex- 
pected to vanish in the chiral limit, whereas the matrix elements {K^\Oi\K^) 
{i = 2, 3, 4, 5) remain finite. For the latter, close to the chiral limit, we expect 
a mild dependence on the quark masses, as given in eq. (3). 

Omitting terms which are of higher order in Chiral perturbation theory 
(xPT), the 5-parameters are usually introduced using the expressions [16] 

3 \m^(/i) +md{n)J 



(ir°|03(/^)|K°) = lf ^^ .X m^^^BM, (5) 

3 ymsifi) +md{lJ,)J 

(ir°|04(/i)|i^°) = 2f /^^ ) MlflB,{^,), 

mOM\K') = I ( . ^^^ . . ) ' M^/i55(/i) , 
3 \ms{fi) +md{iJ,)J 

where Oi{^) and ms{^) + ^^(/i) denote operators and quark masses renor- 
malized at the scale /i. For the four-ferniion operators the most common 
renormahzation schemes are HV and NDR, although DRED is also occa- 
sionally used. 

In (5), the matrix element of the operator Oi is parameterized in terms 
of well-known experimental quantities and -Bi(/i) (Bxifi) = -Bi(//)) has 
been computed with great precision on the lattice [5]- [9]. The expression of 
the matrix elements {K^\Oi\K^) {i = 2,3,4,5) depends, instead, quadrati- 
cally on the quark masses. Therefore, the "conventional" parameterization 
introduces a redundant source of systematic error which can be avoided 
by parameterizing the matrix elements in terms of measured experimental 
quantities 0. 

To overcome this problem, we propose the following new parameterization 
of the AS* = 2 operators 

(ir°|0,(/i)|K°) = Ml.flB,{^i) {i = 2, 3, 4, 5) . (6) 

The advantage of (6) is that it is still possible to work with dimension- 
less quantities, the i?j(/i)s, without any reference to the quark masses. In 
practice, one computes on the lattice the ratios 

~ / ir°|O.M|K°) \ 

\ ^^^K*Jk J latt 

with i = 2,3,4,5, and derive the physical amplitudes, in GeV^, using 

(ir°|0.(/i)|ir°)=5,(Mi./i)^^^. (8) 

In sec. 4 we use the new parameterization to extract the matrix elements 
of the operators Oi {i = 2, 3, 4, 5) from three-point correlation functions. 
Another possible definition is given by the ratios 



^ Most recently, S-parameters with the standard definitions (5) have been com- 
puted in refs. [15,16]. 



from which the physical matrix elements read 



{K'mf,)\K') = -i?.Bi(/i) (M^./^j^^^ • (10) 

In equation (9) the factor M^^/Mj^, has been introduced in order to have a 
smooth behavior of the ratios Ri as a function of the quark masses. 

The Ri may be useful to directly estimate the relative contribution to 
K^-K^ mixing coming from physics beyond the Standard Model [25]-[27]: 

where Ci(/i) and Ci{fi) are the Wilson coefficients of the corresponding op- 
erators. In the Standard Model, the coefficient Ci is the only one different 
from zero. With our data, the error on {K^\Oi{n)\K'^) is rather large. For 
this reason, we only present numerical results obtained with the parameter- 
ization (6). 
A/ = 3/2 operators: 

The study of A/ = 3/2 K ^' rni amplitudes requires the computation 
of the matrix elements (ttttIOj \K) of the following left-left and left-right 
operators [19,20] 

O?^' = s"7^(l - 75)rf" u^-f^il + 75)^'' +{d^u)-{u^d,u^ d) 
08^' = s"7^(l -75)rf^^^7M(l + 75)w° + {d^u)-{u-^d,u^ d){12) 
OT = s"7^(l - 75)^" u^-f.il - 75)w'' +{d^u)-{u^d,u-^d) . 

These matrix elements are important for the calculation of e'/e in the Stan- 
dard Model. In the chiral limit, the (ttttIO/ \K) matrix elements can be 
obtained, using soft pion theorems, from (7r+|0j l-ft'^) (to which only the 
parity-even parts of the operators contribute) 

{nn\0,\K)j=, = --i^(vr+|Or|ir+) . (13) 

The latter can be computed on the lattice using only the three-point cor- 
relation functions. For degenerate quark masses, rris = rrid = m, and in the 
chiral limit, we find 

lim(7r+|0f 17^+) = - lim(^0|O5|ir°) = -M^/^ MmB.ifi) 
\imU+\Of'\K+) = -\im{K'\0,\K') = -Mlfl \imB,{^) 

hm (TT+IOf |ir+) = \Mlfl hm i?,(/.) 

m-^U o m— »U 



where all the matrix elements are parameterized through well known exper- 
imental quantities. Thus we can predict the {'k^\0^ \K^) matrix elements 
from the chiral limit of suitable i?-parameters and from the chiral limit of Bi. 

3 Renormalization Group Invariant Operators 

In this section, we give the main formulae which are necessary to define the 
Renormalization Group Invariant Wilson coefficients and operators for the 
most general A^* = 2 effective Hamiltonian. The procedure is generalizable to 
any effective weak Hamiltonian. 

Physical amplitudes can be written as 

(F|K//|/) = (F|0(^)|/)-(?(/i), (14) 

where Oi^ji) = (Oi(/i), 02(/i), . . . ,OAr(/i)) is the operator basis, e.g. the basis 
defined in (4), and C{fi) the corresponding Wilson coefficients (see for exam- 
ple [18,27]) represented as a column vector. C{fi) is expressed in terms of its 
counter-part, computed at a large scale M, through the renormalization-group 
evolution matrix W[fi, M] 

C{fi) = W[fi,M]C{M). (15) 



The initial conditions for the evolution equations, C{M), are obtained by 
perturbative matching of the full theory, which includes propagating heavy- 
vector bosons (W and Z°), the top quark, SUSY particles, etc., to the effective 
theory where the W, Z^, the top quark and all the heavy particles have been 
integrated out. In general, C{M) depends on the scheme used to define the 
renormalized operators. It is possible to show that W[fi, M] can be written in 
the form 



W[fi,M] = M[fi]U[fi,M]M-^[M] , (16) 

where U is the leading-order evolution matrix 

"«s(M)' 



U[^^,M] 



and the NLO matrix is given by 



tts(/i) 



17g''"/2/3o 

(17) 



M[/x] = i + ^ J[A(/i)] . 



where ■Jq is the leading order anomalous dimension matrix and J[A(/i)] is 
defined in [28] and can be obtained by solving the Renormalization Group 
Equations (RGE) at the next-to-leading order. 

The Wilson coefficients C{^) and the renormalized operators 0(/x) are usually 
defined in a given scheme ( HV, NDR, RI), at a fixed renormalization scale /x, 
and depend on the renormalization scheme and scale. This is a source of con- 
fusion in the literature. Quite often, for example, one finds comparisons of B- 
parameters computed in different schemes. Incidentally, we note that the NDR 
scheme used in the lattice calculation of ref. [15] differs from the standard NDR 
scheme of refs. [19,20]; on the other hand, the HV scheme of refs. [19] is not 
the same as the HV scheme of refs. [20] . In some cases, the differences between 
different schemes may be numerically large, e.g. -Bg ~ 1.3 -Bg • -^^^ 

these reasons, the standard procedure is not entirely satisfactory, especially 
when the (perturbative) coefficients and the (non-perturbative) matrix ele- 
ments are computed using different techniques, regularization, schemes and 
renormalization scales. To avoid all these problems, we propose a Renormal- 
ization Group Invariant (RGI) definition of Wilson coefficients and composite 
operators which generalizes what is usually done for Bk, by introducing the 
RGI 5-parameter B^ and for the quark masses [29,30]. The procedure is 
straightforward: from eq. (17), we define 



II I 71 4- I II / \ I ' 

w 



M[f,] [«,(/i)]-^O^"/2/30 ^ (19) 



which, using eqs. (16) and (19), gives 

W[fi,M]=w~^[l2]w[M] . (20) 

The effective Hamiltonian (14) can then be written as 

Ueff = 6{i^) ■ C{fi) = 6{fi)W[i2, M]C{M) 

= d(/i)t&-^[/i] ■ w[M]C{M) = O^^^ ■ (7^^^ , (21) 

with 

C^^^ = w[M]C{M) 

d^G^ = d(/i)-t&~^[/i] . (22) 

(jRGi ^j^^ qRGI g^j,g scheme and scale independent at the order at which the 
Wilson coefficients have been computed (NLO in most of the cases, NNLO for 
quark masses). 



The S-paraineters defined in eqs. (6) satisfy tlie same renormalization group 
equations as tfie corresponding operators. Tlie RGI 5-parameters are tlien 
obtained from the relation 

B^ = T.BA^^)w{^^)J.\ (23) 



4 Numerical results 

In this section we present numerical results for the matrix elements of the 
basis (4) and for the electro-penguin operators (12), obtained using our new 
parameterization. 

All details concerning the non-perturbative renormalization of lattice oper- 
ators and the extraction of matrix elements from correlation functions have 
been presented elsewhere and are not repeated here. The reader can find them, 
for example, in refs. [14,16], (see also ref. [31] for a complete discussion of the 
non-perturbative renormalization techniques for AF = 2 operators and for 
references). In particular in [16], the same set of lattice data was analyzed and 
is compared here with the results of the our new study. 

Let us start by giving a simple argument which shows the correlation existing 
between the values of the "conventional" 5-parameters and the quark masses. 
In order to extract the operator matrix elements the following two- and three- 
point correlation functions are used 

Gpit^^p) = ^(P(a;)pt(o))e-P"^^ , ^^(ix,^) = E(^o(a:)Pt(0))e-P^ 

X X 

Goitx, ty;p, q) = ^(pt(j/)O(0)pt(a;))e-P-^V^ (24) 

where x = (x, tx),y = {y, ty) and P, Aq and O stand for the renormalized pseu- 
doscalar density, the fourth component of the axial current and four-fermion 
operator respectively 0. By forming suitable ratios of the above correlations 
and looking at their asymptotic behavior at large time separations, one can di- 
rectly obtain the P-parameters defined in eqs. (5). At the leading order in the 
chiral expansion, or using the definitions of ref. [16], one finds (for z = 2, 3, 4, 5) 



«. . ,5%%fL ^ »™> ^ con... X B.(., . ,25, 

Gp{U,p)Gp{ty,q) |(0|P|is:")|^ 

On the other hand, quark masses are extracted using vector or axial-vector 
Ward identities, e.g. 



^ In this study, all correlation functions are evaluated with degenerate quark 
masses. 



Table 1 

Matrix elements at the renormalization scale // = a^^ ~ 2 GeV, corresponding to 
fi'^a'^ = 0.96 and //^a^ = 0.62 at f5 = 6.0 and 6.2 respectively. All results are in 
the RI (MOM) scheme. In the first two columns the results of the present study 
obtained with the new parameterization are given. In the last two columns we show 
the results obtained with the "conventional" method in ref. [16] on the same set of 
data (for (Oi) the two parameterizations coincide). (O^g ) stand for (vr"^|07 8 \K'^) 

3/2 

with Oj g given in eq. (12). 





New 


Old 


m 

fi ~ 2GeV 


/3 = 6.0 
this work 


/3 = 6.2 
this work 


/3 = 6.0 
[16] 


/3 = 6.2 
[16] 


(Oi) 
Bi 


0.012(2) 
0.70(15) 


0.011(3) 
0.68(21) 


0.012(2) 
0.70(15) 


0.011(3) 
0.68(21) 


(O2) 
B2 


-0.079(10) 
0.72(9) 


-0.074(8) 
0.67(7) 


-0.073(15) 
0.66(3) 


-0.073(15) 
0.66(4) 


(O3) 
Bs 


0.027(2) 
1.21(10) 


0.021(3) 
0.95(15) 


0.025(5) 
1.12(7) 


0.022(5) 
0.98(12) 


(O4) 
Bi 


0.151(7) 
1.15(5) 


0.133(12) 
1.00(9) 


0.139(28) 
1.05(3) 


0.133(28) 
1.01(6) 


(O5) 
B5 


0.039(3) 

0.88(6) 


0.029(5) 
0.66(11) 


0.035(7) 
0.79(6) 


0.029(7) 
0.67(10) 


{oT) 


-0.019(2) 
0.65(5) 


-0.011(3) 
0.38(11) 


-0.020(5) 
0.68(7) 


-0.014(5) 
0.46(13) 


{oT) 

bT 


-0.082(4) 
0.92(5) 


-0.068(8) 
0.77(9) 


-0.092(19) 
1.04(4) 


-0.087(19) 
0.98(8) 



ms{ij) +md{it) 



{0\d,A,\K'>) 
{0\P\KO) 



(26) 



where the quark masses ms{^)+md{ii) and the renormalized pseudoscalar den- 
sity P are given, by definition, in the same renormalization scheme. This hap- 
pens because matrix elements of the "good" axial currents, and consequently 
the product {ms{fi)+m(i{ii))P are regularization, renormalization scheme and 
scale independent. In practice, for a given renormalization scheme and scale 
adopted to renormalize P, the matrix elements {0\d^A^\K^) and {0\P\K^) are 
computed numerically and from their ratio one obtains ms{ij) + md{ij). 
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Fig. 1. (^^1041^^") in GeV^ computed with the new parameterization as a func- 
tion of the sum of the quark masses rng + rud renormalized in the MS scheme at 
fi =~ 2 GeV. The quark masses span a range of values compatible with theoretical 
estimates. 



0.2 



0.15 




0.05 



150 



m„ 



m. 



From a comparison of eqs. (25) and (26), it is evident that, for large values 
of (0|P|i^''), small values of the quark masses and of the 5-parameters will 
be simultaneously obtained. The ratio Bi/{(ms{fi) + md{fi))'^ , corresponding 
to the "physical" matrix elements, including the mass factors appearing in 
eqs. (5), will however be much less dependent on the quark masses since Bi 
and ms{[j) + md{[j) are strongly correlated, i.e. (0|P|i^°) oc l/rriq. It would 
be necessary that calculations of "conventional" 5-parameters provide at the 
same time the value of the quark masses obtained in the same calculation. This 
is equivalent to give, as suggested in eq. (6), the matrix elements in physical 
units. 



In this work we have used the lattice correlation functions computed in [16]. 
They have been obtained from numerical simulations at /? = 6.0 (460 con- 
figurations. Run A) and 6.2 (200 configurations. Run B) with the tree-level 
Clover action, for several values of the quark masses and for different me- 
son momenta. We have used the "lattice dispersion relation" sinh [—2) ~ 

sinh^ ( T ) + S j=i,3 sin^ ( ^ ) ■ The physical volume is approximatively the same 
on the two lattices. Statistical errors have been estimated with the jackknife 
method. 
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The main numerical results that we have obtained are listed in table 1. They 
have been computed in the RI(MOM) scheme at /i = 2 GeV as follows: 

• In the first two columns of the table the results, obtained with the new 
parameterization, are given. The 5-parameters are obtained by extrapolat- 
ing/interpolating the ratios 

^ Ml GoXt.,ty,p,q) {K\q)mK\p)) 

' M^GA{t.,p}GA{ty,q) M^fl ■ ^ ^ 

to the physical point using the lattice-plane method [32]. Mp and My are 
the masses of the pseudoscalar and vector mesons computed at the same 
quark masses as the correlation functions. Then the matrix elements (Oj) in 
GeV^ are obtained from eq. (8). The standard BiS given in these columns 
are obtained from eqs. (5) by using the matrix elements in GeV^ in the same 
column and a "conventional" quark mass {ins + 'm(i)'^^{'2, GeV) = 130 MeV. 
This is the value of the strange quark mass that we have obtained on the 
same sets of data [33] . 

• In the last two columns we show the results obtained with the "conven- 
tional" method in ref. [16] on the same set of data. In this case we first ob- 
tain the 5-parameters, defined as in eq. (5), using the procedure described 
in [16]. The matrix elements in GeV"^ reported in these columns are derived 
from the "conventional" B-parameters in eqs. (5) with the same quark as 
before of (m^ + md)^^{2 GeV) = (130 ± 15) MeV. In this case we have to 
include the uncertainty coming from the error in the determination of the 
quark mass. The errors on the 5-parameters and the quark masses are con- 
sidered as independent to mimic what is usually done in phenomenological 
analyses. 

The results in table 1 show the convenience of the approach proposed in this 
paper. Most phenomenological analyses, which use the standard procedure, 
adopt values of 5-parameters and quark masses taken from different determi- 
nations. For example the 5-parameter is taken from a calculation at a fixed 
value of the lattice spacing in the quenched approximation whereas the quark 
mass is taken from some compilation of lattice results extrapolated to the con- 
tinuum and including some evaluation of the quenching errors. In other cases, 
the value of the 5-parameter is a sort of average biased by results from the 
large N approach and lattice calculations and the mass is taken from an av- 
erage from lattice and QCD sum rules. In this way, one misses the correlation 
between the value of the S-parameter and the mass. As a result, the uncer- 
tainty on the physical matrix element is much larger. In figure 1 we present 
the values of the matrix element {K^\0^\K^) obtained with the new param- 
eterization for different choices of m^ -|- rud, in a range of values compatible 
with theoretical estimates. To be more specific we have computed {K^\04\K^) 
for different values of rUs + nid, where m^ and m^ are the masses of the quarks 
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in the meson. This means that different choices of rris + rrid also correspond to 
different values of Mj^o (and of fx)- 

Although we have data at two different values of the lattice spacing, the sta- 
tistical errors, and the uncertainties in the extraction of the matrix elements, 
are too large to enable any extrapolation to the continuum limit a ^ : 
within the precision of our results we cannot detect the dependence of B- 
parameters on a. For this reason, we estimate the central values by averaging 
the 5-parameters obtained with the physical mass M^^ at the two values of 
f3. Since the results at /5 = 6.0 have smaller statistical errors but suffer from 
larger discretization effects, we do not weight the averages with the quoted 
statistical errors but simply take the sum of the two values divided by two. 
As far as the errors are concerned we take the largest of the two statistical 
errors. This is a rather conservative way of estimating the errors. In order to 
compare the results of Run A and Run B, we have chosen the same physical 
renormalization scale /i. Using estimates of the lattice spacing (a^^ = 2.12(6) 
at /3 = 6.0 and a^^ = 2.7(1) at /5 = 6.2) of ref. [33], we have taken fi^a^ = 0.96 
and fi'^a^ = 0.62, corresponding to fj, = 2.08 GeV and /i = 2.12 GeV, at 
P = 6.0 and 6.2 respectively. We quote the results as obtained at /x = 2 GeV, 
since the running of the matrix elements between /i ~ 2.1 and 2.0 is totally 
negligible in comparison with the final errors. 

5 Physical Results 

Our best estimates of the matrix elements of AS* = 2 operators in the RI 
scheme at /x = 2 GeV are given in table 2. Note that these matrix elements 
are enhanced by a factor ~ 2 -^ 12 with respect to the SM one ((Oi)). For this 
reason, K^ — K^ mixing is a promising observable to detect signals of new 
physics at low energy [18,27]. 

From the matrix elements obtained non-perturbatively in the RI scheme, we 
derive the RGI 5-parameters using continuum perturbation theory at NLO. 
These have been computed from eqs. (23) with Q;s(2GeV) = 0.31, correspond- 
ing to as{Mz) = 0.118, evaluated with the appropriate number of active 
flavors ( rij = 4 at 2 GeV). The results are given in table 2. This choice mim- 
ics what is usually done in phenomenological analyses which use lattice QCD 
estimates. It corresponds to the assumption that the results in the Rl-scheme 
are the "physical" matrix elements, up to some undetermined unquenching 
errors. 

Our best estimate of the matrix elements of A/ = 3/2 operators in the RI 
scheme at /i = 2 GeV are 

{ti+\0'^^'^\K+) = -(0.015 ± 0.004)GeV^ , 

{7r+\0f^\K+) = -(0.075 ± 0.008)GeV^ . (28) 

Since the analyses of refs. [24], [34]- [39] are done using Wilson coefficients in 
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Table 2 

Matrix elements in GeV^ at fi = 2 GeV in the RI scheme and their RGI values with 
-/=4 



a 



m 


RI 


RGI 


(Oi) 


0.012(3) 


0.017(4) 


(02) 


-0.077(10) 


-0.050(7) 


(03) 


0.024(3) 


0.001(7) 


(04) 


0.142(12) 


0.068(6) 


(05) 


0.034(5) 


0.038(5) 



Table 3 

Matrix elements in GeV^ at fi 



2 GeV in the NDR and HV schem,es and with 



as 



{2GeV) = 0.31. 



{oT) 


NDR 


HV 


(vr+iof |i^+) 


-0.021(4) 


-0.033(5) 


{7it:\0-j\K)i=2 


0.11(2) 


0.18(3) 


{^^\0T\K+) 


-0.095(10) 


-0.114(12) 


{Tn:\Os\K)i=2 


0.51(5) 


0.62(6) 



NDR and HV, by using the matching coefficients between the RI and these 
schemes [24] 



{Oi 



3/2\ 



NDR.HV 



Jij 



^MAO A NDR,HV 
47r ^'^ 



(Of) 



RI 



(29) 



where 



Ar, 



NDR 



| + iln2 -2 



2 In 2 



Ar 



HV 



-8-21n2 



2-21n2 -f + |ln2 



-| + iln2 



-2-21n2 -f + |ln2 



and the results in (28), we have computed the matrix elements (tt^IO/ \K^) 
given in table 3. To obtain (7r7r|Oj|i^)/=2 we have used the chiral relation of 
eq. (13). This entails further uncertainty in the numerical evaluation of the 
physical matrix elements: within our accuracy, we may use in eq. (13) fx 
instead of f^^. Moreover, to obtain the matrix elements in the MS scheme 
from those obtained non-perturbatively in the Rl-scheme, we have chosen the 
unquenched value Q;s(2GeV) = 0.31 but, within the quenched approximation, 
we could have chosen the quenched value of a^ as well. We estimate that, due 
to these effects, the final error is about twice that quoted in table 3. 
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It is interesting to compare our result for {tx'^\0-j \K^) in the NDR scheme 
with the recent estimate from large A^ expansion [40] . The two determinations 
differ by more than two a with respect to the error quoted in table 3. 

6 Conclusions 

In this work we have introduced a new parameterization of four fermion oper- 
ator matrix elements which does not involve quark masses and thus allows a 
reduction of systematic uncertainties. As a result the apparent quadratic de- 
pendence of e'/e on the strange quark mass is removed. We have also defined 
Renormalization Group Invariant matrix elements to simplify the matching 
between the lattice and continuum renormalization schemes. We have used 
these definitions to compute matrix elements of A/ = 3/2 and SUSY AS = 2 
four fermion operators on the lattice in the quenched approximation. The sim- 
ulations have been performed at two different values of the lattice spacing and 
the renormalization constants of the operators are calculated non perturba- 
tively. 
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